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Abstract 

We study experimentally the shape of a granular heap formed pouring 
a granular material into a vertical Hele-Shaw cell and analyze the effect of 
the grain impact energy. We propose a continuous model for the steady 
profile of the heap that explicitly considers energy dissipation of flowing 
grains through inelastic collisions. We solve the model analytically and 
analyze the resulting height profile as a function of several parameters, 
such as the restitution coefficient of the grains and their impact energies. 
We find good agreement between theory and experiments. 



1 Introduction 

The study of granular materials has received wide attention in the last few years 
[|l| ||, |!| , because of the challenging problems posed by this state of matter. Sand 
heaps are among the most familiar granular systems that have been studied ex- 
perimentally and theoretically. Despite the apparent simplicity a sand heap 
presents a number of intriguing phenomena, such as avalanches Q internal tex- 
ture j|] , segregation and stratification |(| . Several theoretical models have been 
employed to describe the dynamics and the static profile of a sandpile, ranging 
from purely numerical models, such as lattice gas [Q |^] and molecular dynamics 
simulations || , to discrete lattice models [[[o) and continuous equations (ll], [l2| . 

Granular heaps can be generated by slowly pouring grains (i.e. spherical 
beads, sand) from a given height over a rigid support (3D piles) or into a ver- 
tical Hele-Shaw cell (2D piles). The pile typically displays a well defined slope, 
corresponding to the angle of repose, with clear deviations at the top and at the 
tail. In particular, the tail of a 2D pile was found experimentally to follow a 
logarithmic dependence |u| This result was anticipated by the solution of 
a continuous model solved in a 2D silo geometry jl2| and by theoretical argu- 
ments jl3| . The angle of repose depends on the characteristics of the material 
such as density, humidity, packing history and boundary conditions. 
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Here we analyze the shape of the top of the granular heap, focusing on the 
effect of the impact energy of the particles. To this end, we perform a set of 
experiments with glass and leads beads poured at very slow rate into a vertical 
Hele-Shaw cell and study the changes of the heap profile when grains are poured 
from different heights. We propose a simple continuous model to describe the 
steady-state properties of the grains flowing on the top of the pile, providing 
an explicit expression for the heap profile. The model treats explicitly the 
dissipation of the kinetic energy of the flowing grains due to inelastic collisions 
and provides a good fit to the experimental data on glass and lead beads. 

2 Experiments 

Experiments are performed within a vertical Hele Shaw cell (20x30cm) with a 
fixed thickness of 5mm Jig] . Lead beads, with diameter d — 0.2cm are poured 
through a funnel placed above the center of the cell. We carefully determine 
the impact energy of the particles , controlling the distance Zq between the top 
of the heap and the funnel: the impact energy is thus e = mgZ . The falling 
heights range from Zq = 1cm to Zq = 14cm and are determined with a with 
a precision of AZ = 0.5cm. Experiments are performed at very low particle 
flux, dropping one particle at a time through the funnel. In the range of falling 
heights employed [|l6| , the reorganization of the pile after impact is usually very 
weak and the particle jumps down the slope and is eventually trapped, giving 
rise to a local change in the slope of the pile. The dynamics of the beads can 
be sometimes very complex: we observe collisions with the walls of the cell, 
rolling motion, ejection of other particles and occasionally big avalanches of a 
few layers. 

After a short transient stage, the heap grows steadily and the profile trans- 
lates at constant velocity, proportional to the input flux. In this regime we 
retrieve the experimental heap profiles (see Fig. |l|)from image analysis and we 
average them. The averaging is done superimposing the profiles obtained at 
different time steps in the steady state. In Fig. || we show the averaged pro- 
files for different impact energies. It can be observed that the heap has a flat 
top and shows a well defined angle of repose, which turns out to be essentially 
independent on Zq, at least for the experimental values used. When the tails 
of the heaps obtained with different Zq are superimposed, we observe that the 
maximum height h(0) decreases linearly with Zq (inset of Fig. ||). We find that 
the shape of the heap is well described by the form 

h{x)=A-6 c x-6 c &- x /t\ (1) 

where 9 C is the angle of repose, £ is a correlation length and A is an inessential 
constant, that can be eliminated through an appropriate vertical translation of 
the profile. From Eq. ([!]) we have that h(0) = A — 9 C ^, which implies that £ is 
proportional to Zq, which is verified by the fitting parameters (see Fig. ||) This 
result reflects the fact that Zq is, together with the particle diameter, the only 
characteristic length present in the system. 
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The basic features of these results are quite robust and do not depend on the 
particular material used. We verified this by forming heaps with glass beads of 
considerably smaller diameter [d = 0.025cm), poured at low flux W = Q.32g/s. 
In this case Zq can be controlled very precisely by moving the outlet up at the 
same velocity as the growth of the top of the heap. In this experiment, grains 
are poured from one of the two sides of the Hele-Shaw cell. We see in Fig. [|that 
also in this case h(0) is linear in Zq and the decay is reasonably well fitted by an 
exponential, excluding from the fitting region the first 0.5cm, which correspond 
to the radius of the funnel. Finally in this case, the measured angle of repose 
displays a small dependence on Z , of the order of A8 = 0.05. Looking at Fig. ||, 
we note that the linear dependence of £ from Z is steeper in the case of glass 
beads by roughly a factor 2. The origin of this effect is probably due to the 
different restitution coefficients of the materials. 



In order to understand the experimental observations, we consider a sandpilc 
model in which energy dissipation is explicitly treated. The complete dynamical 
evolution of all the grains and their energy dissipation is extremely difficult to 
treat explicitly. The evolution of the heap is indeed a very complex phenomenon, 
being due to the way momentum is transfered during particle collisions, fric- 
tional properties and other forms of energy dissipation. However, we note that 
in our experiment the growth of the heap displays a remarkable stationary state, 
which we can try to describe with a simplified model taking into account essen- 
tially only energy dissipation. 

Alonso et al. |lC[ ] have introduced a simple automaton model which has been 
used to derive the angle of repose of a pile in terms of microscopic parameters, 
such as the restitution coefficient of the grains. In this model, grains are poured 
one by one on the top of a heap with an initial kinetic energy eo. As the grains 
fall down the pile, parameterized by the local height h(x), their kinetic energy 
e decreases due to inelastic collisions and increases because of gravity as 



where m is the mass of the grains, 5 is the characteristic size of the jumps, g is 
the gravity acceleration and r is an "effective" restitution coefficient. The grains 
jump to the right until their kinetic energy is less than a threshold u. When a 
grain comes to rest at x the height h of the column x is increased by one unit 
and a new grain is poured at the origin. This model has been simulated and, 
despite its simplicity, was shown to explain several experimental results Jl7j . 

We study here a continuous version of the model in which the steady-state 
kinetic energy of the rolling grains follows Eq. (0) 



steady-state energy profile in terms of a given slope profile 7(2;). To close the 



3 Model 



e(x + ft) = r(e(x) + mg(h(x) — h(x + 5))), 



(2) 



where 7(2;) 
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problem, we need to specify how the steady-state slope profile depends on the 
energy profile. We assume that the system in the steady state is at the verge 
of stability and that the variations of the local slope compensate the difference 
between the kinetic energy e and the trapping energy u: the slope will decrease 
for e > u and increase for e < u. The simplest equation of this form is given by 

S^- = T(u-e(x)), (4) 

where F is a phenomenological parameter. Eq. (Q) implies that for large x the 
grains are at the threshold of mobility e(x) ~ u. 

Finally, we impose the boundary conditions e(0) = cq = mgZo and 7(0) = 0. 
We note that by dividing each term of the equation by mg all the quantities 
can be expressed in terms of lengths. Eqs. (|^-(|) are linear and can be solved 
in order to obtain an explicit expression for the profile of the pile. It is con- 
venient to transform the pair of equations into an homogeneous one using the 
transformation 

E(x) = e(x) - u, G(x) = 7(2) + (I - r)u/r. (5) 
The solution of the homogeneous equations 

S ^M. = (r-i)E{x)+rG(x), (6) 
dx 

&® - -rsw. (7) 



is given by 



e(x) = u + E + e x+x/s + E^e x - x/s (8) 
j(x) = (r - l)u/r + G+e x + x/s + G-e x - x/s (9) 



where A± = (r — 1 ± yjr — l) 2 — 4rr)/2 are the eigenvalues of the matrix 

c -( r -r;) < 10 » 

The coefficients E± and G± can be explicitly obtained imposing the boundary 
conditions. 

In the limit of large x the slope of the pile is given by 700 = (r — l)u/r, from 
which we obtain the angle of repose 6 C = (1 — r)u/ (mgSr). Since the real part of 
A± is negative for all the values of the parameters the model has a well defined 
angle of repose. The relaxation towards the angle of repose is expressed by the 
sum of two exponentials with characteristic lengths £± = 5/X±. Since £ + > 
the relaxation is dominated by £ + that we can identify with £ of Eq. |l|. Finally, 
integrating the equation for 7(2:), we obtain an analytical form for the shape of 
the pile that can be compared with experiments 

h(x) — h(0) — (1 — r)ux/(mgSr) — 
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5[G + /X+(e x + x ^ - 1) + G-\-(e x - x/s - 1)]. 



(11) 



The exponential relaxation to the angle of repose is experimentally observed, 
although in the present form the model can not account for the impact energy 
dependence of the relaxation length. 

An original assumption of Eq. (||) was that the jump length 5 is constant 
during the grain motion. The problem of a ball rolling and jumping over a 
rough inclined plane has received a large attention in recent years |l8|, ^9|. It 
has been experimentally observed that as a function of the angle of the inclined 
plane the motion is decelerated or accelerated, with a region in between where 
the motion is effectively viscous. The energy relaxation was recently studied in 
Ref. |is|| , where the stopping length in the decelerated regime was measured as 
a function of the initial energy. The result clearly indicates that the stopping 
length is linearly dependent on the initial energy. We can show this fact by a 
simple argument. Consider an inelastic ball hitting a plane with energy eo and 
let a be the angle between the plane and the particle velocity before the impact. 
We consider for simplicity that normal and tangential restitution coefficients are 
both equal to r so that at each jump n the energy decreases by e„ = re n -\. 
A straightforward calculation shows that the jump length decreases with n as 
S n = 2e r™ sin(2a)/(mgf), which yields e(x) = e [l — mgx(l — r)/(2e sin(2a))] = 
e (l — x/£). We see here that the characteristic length for energy relaxation £ 
is linear in eo, in agreement with the experimental result seen in Fig. ^| 

4 Conclusions 

In conclusion, we have studied the effect of impact energies on the shape of the 
sand heap. Experiments show an exponential relaxation of the slope to the angle 
of repose, with a characteristic length that is linearly dependent on the impact 
energy of the grains. We model this process by simple relaxation equations and 
explain theoretically the experimental results. 

S. Z. is supported by EC TMR Research Network contract ERBFMRXCT960062. 
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Figure 1: An example of a heap obtained experimentally with lead beads. 
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Figure 2: The height profiles of the heaps obtained experimentally with lead 
beads for different impact energies Zq fitted with Eq. [l| The dashed line has a 
slope 8 C ~ 0.54. In the inset we show the linear behavior of h(Q) with Zq. 
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Figure 3: Linear relation between the characteristic lenght £ and the impact 
energy Zq for lead and glass beads. 
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Figure 4: The height profile obtained in the experiment with glass beads, for 
different impact energies Zq at the lowest available input flux W = Q.32g/s. 
The lines are the fit with Eq. [j]. In the inset we show the linear dependence of 
h(0) with Zq. 
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